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1.  INTRODUCTION  AND  SUMMARY 


1.1  OBJECTIVE 

The  objective  of  this  program  is  to  determine  what  influences  the  sensi¬ 
tivity  of  optical  fiber  to  radiation  and  use  this  knowledge  to  develop  fiber 
with  improved  radiation  resistance.  The  effort  is  directed  toward  graded  in¬ 
dex  fibers  which  are  required  for  high  bandwidth  applications.  Because  of 
the  dopants  in  the  core  which  are  required  to  grade  the  refractive  index, 
these  fibers  are  much  more  sensitive  to  radiation-induced  loss  than  the  lower 
bandwidth  pure-silica-core  fibers. 

The  program  is  based  upon  the  well  established  concept  that  radiation- 
induced  absorption  is  governed  by  pre-existing  defects  in  the  glass.  The  in¬ 
duced  loss  is  due  to  color  centers  which  form  when  existing  defects  trap  the 
electrons  and  holes,  generated  by  ionizing  radiation.  The  emphasis  is  on  ob¬ 
taining  a  fundamental  understanding  of  the  relationship  between  defects  in  the 
as-drawn  fiber  and  the  glass  composition  and  processing  conditio  s.  Using 
this  knowledge,  fibers  with  parameters  predicted  to  be  favorable  will  be  fa¬ 
bricated  and  tested.  For  the  most  radiation -resistant  optical  fibers  devel¬ 
oped,  complete  documentation  on  compositions  and  fabrication  processes  will  be 
provided  in  sufficient  detail  to  allow  transfer  to  a  manufacturing  unit. 


1.2  APPROACH 


Induced  absorption  occurs  when  existing  defects  in  the  glass  form  color 
centers  by  trapping  electrons  and  holes  generated  by  ionizing  radiation. 
Three  technical  approaches  are  being  pursued  which  are  based  upon  the  prin¬ 
ciple  of  reducing  the  concentrations  of  these  defects  or  mitigating  their  ef¬ 
fect. 

I .  Defect  Passivation 

One  method  to  prevent  defects  in  as-drawn  fiber  from  capturing  carriers 
and  forming  color  centers  is  to  transform  them  into  benign  defects.  The  lat¬ 
ter  are  defined  as  either  defects  which  form  color  centers  that  absorb  out¬ 
side  the  wavelength  range  of  interest,  or,  preferably,  ones  without  the  abil¬ 
ity  to  capture  an  electron  or  hole. 

II.  Defect  Minimization 

Reducing  the  concentration  of  defects  in  the  as-drawn  fiber  will  reduce 
the  amount  of  color  center  absorption  following  irradiation.  Common  sources 
of  defects  are  dopants,  mechanical  stress,  and  thermomechanical  stress. 

III.  -  Recombination  Centers 

Many  defects  and  impurities  have  localized  electronic  states  located  deep 
within  the  energy  gap  of  the  glass.  These  will  serve  as  recombination  cen¬ 
ters  if  the  states  have  the  ability  to  trap  both  electrons  and  holes.  These 
centers  will  reduce  induced  loss  if  they  can  compete  successfully  with  other 
defects  for  the  carriers  generated  by  the  radiation. 

The  program  is  divided  into  three  phases.  In  the  first  or  Exploratory 
Phase  all  of  the  three  technical  approaches  will  be  pursued  simultaneously. 
This  phase  has  a  duration  of  *15  months.  The  second  or  Refinement  Phase 
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will  last  >10  months  during  which  the  activity  will  focus  on  only  the  most 
promising  approach  or  combination  of  approaches.  The  final  or  Optimization 
Phase  involves  fine  tuning  of  the  successful  approaches,  including  the  explo¬ 
ration  of  tradeoffs.  This  report  covers  the  first  12  months  of  the  Explora¬ 
tory  Phase. 


1.3  PROCEDURE 

The  program  involves  both  an  experimental  and  a  theoretical  effort.  The 
former  consists  of  carefully  chosen  experiments  based  upon  established  con¬ 
cepts  from  the  physics  of  defects  in  glasses.  The  primary  experimental  tool 
is  the  measurement  of  radiation-induced  loss  in  the  fibers.  This  is  being 
done  both  In  situ  at  RADC  at  a  fixed  wavelength,  and  as  post  irradiation 
spectral  measurements  at  GTE  Laboratories.  Differential  modal  attenuation  is 
also  being  applied  in  some  of  the  approaches.  More  selective  techniques  such 
as  Raman  scattering  and  photoluminescence  are  also  being  tried,  and  the  lat¬ 
ter  method  shows  great  promise  for  characterizing  defects  in  fibers.  These 
experimental  methods  are  being  supplemented  with  molecular  dynamic  simula¬ 
tions  of  the  glass  to  determine  defect  types  and  concentrations.  This  in¬ 
volves  computer  generation  of  glass  structures  which  are  then  analyzed  to 
identify  defects.  The  objective  of  the  molecular  dynamics  simulations  is  to  es¬ 
tablish  the  correlation  between  structural  defects  and  the  dopants  used  either 
to  modify  the  index  of  refraction  or  to  introduce  recombination  centers. 


1.4  SUMMARY  OF  RESULTS 


All  of  the  approaches  investigated  to  date  have  been  extremely  effective 
at  reducing  induced  loss.  Fibers  have  already  been  produced  which  are  com¬ 
parable  in  radiation  hardness  to  the  best  doped-core  fibers  reported  in  the 


literature.  No  results  are  yet  available  for  the  recombination  center 
approach,  but  combinations  of  at  least  the  defect  passivation  and  defect  mini¬ 
mization  approaches  now  appear  to  be  extremely  promising  avenues  to  further 
improvements. 

In  the  defect  passivation  approach,  the  major  effort  has  been  in  the  area 
of  hydrogen  treatment  of  fibers  to  convert  existing  defects  into  types  which 
result  in  less  attenuation  when  irradiated.  This  technique  is  showing  great 
promise  for  hardening  fibers.  The  first  series  of  experiments  were  conducted 
on  fibers  with  a  Ce/P-codoped  core.  The  fibers  which  had  been  treated  with 
hydrogen  prior  to  irradiation  demonstrated  significantly  less  induced  loss. 
The  effect  was  found  to  depend  upon  the  extent  of  the  reaction  of  hydrogen 
with  the  defects:  treatment  with  a  higher  hydrogen  partial  pressure  im¬ 
proved  the  hardness  still  further.  It  was  also  found  that  hydrogen  treatment 
following  irradiation  significantly  reduced  the  induced  loss  and  made  the  fi¬ 
bers  much  harder  in  subsequent  irradiations.  Since  phosphorus -doped  fibers 
are  very  radiation  sensitive,  this  technique  will  now  be  applied  to  much 
harder  Ge/F-doped  fibers  discussed  below. 

An  effort  was  made  to  find  a  characterization  technique  with  greater  se¬ 
lectivity  than  attenuation  measurements  for  the  identification  of  specific  cen¬ 
ters.  Raman  scattering  and  photoluminescence  investigations  of  the  fibers 
were  undertaken,  and  although  Raman  scattering  does  not  seem  to  be  sensi¬ 
tive  to  the  centers  of  interest,  photoluminescence  was  found  to  be  an  effec¬ 
tive  tool.  Significant  differences  were  found  among  fibers  subjected  to  differ¬ 
ing  treatments.  In  particular,  the  intensity  of  the  luminescence  band  at 
6S0  nm,  attributed  to  drawing  and/or  radiation  induced  centers,  was  found  to 
have  an  inverse  correlation  with  loss.  The  luminescence  intensity  also  in¬ 
creased  with  photobleaching.  These  results  support  our  premise  that  there 
are  a  wide  variety  of  defects  present  in  fiber,  some  more  deleterious  than 
others  in  a  radiation  environment,  and  that  defects  can  be  converted  from  one 
type  to  another. 


The  defect  minimization  approach  was  the  other  area  which  was  focused 
on  in  the  current  reporting  period.  The  molecular  dynamic  simulations  were 
successful  in  generating  pure  silica  glass,  pure  germania  glass,  and  a 
10  mole%  Ge-doped  silica  glass.  The  glass  simulation  is  now  well  in  hand  and 
the  emphasis  has  turned  to  the  analysis  programs,  particularly  the  one  which 
calculates  the  radial  distribution  functions  and  coordination  numbers.  This 
analysis  is  now  being  applied  to  the  glasses  to  determine  the  types  and  con¬ 
centrations  of  structural  defects.  The  Ge-doped  silica  has  been  the  most 
thoroughly  analyzed  glass,  and  no  evidence  was  found  for  defects  resulting 
from  anomalous  silicon-oxygen  distances.  More  detailed  analysis  of  these 
glasses  will  be  required. 

Important  results  have  been  obtained  in  the  investigation  of  the  role  of 
drawing-induced  defects  in  radiation  sensitivity.  Four  fibers  from  the  same 
preform  were  drawn  under  different  conditions  of  temperature,  speed,  and 
tension.  Clear  evidence  was  found  that  a  lower  draw  temperature  produces 
harder  fiber.  No  unambiguous  correlation  could  be  found  for  the  effect  of 
drawing  speed.  These  fibers  had  a  Ge/F-codoped  core  and  a  F-doped  de¬ 
pressed  cladding.  The  induced  loss  was  quite  low  for  all  these  fibers  and 
was  comparable  to  the  best  values  reported  for  doped-core  multimode  fiber  in 
the  literature.  The  conditions  for  low  drawing-induced  defect  concentrations 
established  by  this  study  are  being  applied  in  the  other  areas  of  the  pro¬ 
gram. 

The  investigation  of  core/clad  interfacial  stress  used  differential  modal 
analysis  in  addition  to  standard  loss  spectrum  measurements.  The  DMA  meas¬ 
urements  revealed  that  prior  to  irradiation  the  higher  modes  suffered  greater 
attenuation  in  the  fibers  with  the  highly  stressed  core/clad  interface  than  in 
the  low  stress  fibers.  After  irradiation  the  differences  between  the  high  and 
low  stress  fibers  were  found  to  be  insignificant  from  both  attenuation  and 
DMA  measurements.  It  was  concluded  that  core/clad  interfacial  stress  does 
not  play  a  major  role  in  determining  the  radiation  sensitivity  of  a  fiber  since 
only  the  higher  modes  are  affected  and  these  modes  do  not  propagate  long 
distances  under  any  circumstances. 


Work  has  recently  begun  on  improving  radiation  resistance  by  reducing 
the  dopant  concentration  in  the  core.  Since  doping  introduces  defects  into 
glass  by  altering  the  structure,  we  expect  that  induced  loss  can  be  reduced 
by  lowering  the  dopant  concentration  or  by  using  combinations  of  dopants 
which  do  not  introduce  defects.  To  this  end,  three  preforms  have  been  fa¬ 
bricated  with  different  concentrations  of  germanium  and  fluorine  in  the  core. 
The  fibers  have  been  drawn  and  are  in  the  process  of  being  tested. 


2.  DEFECT  PASSIVATION 

2.7  HYDROGEN  TREATMENT 

One  method  to  prevent  defects  in  as-drawn  fiber  from  capturing  carriers 
and  forming  color  centers  is  to  transform  them  into  benign  defects  which  do 
not  absorb  in  the  spectral  region  of  interest.  This  approach  has  already 
been  exploited  in  reducing  radiation  damage  in  bulk  vitreous  silica  through 
the  use  of  hydrothermal  treatment  [1].  The  effect  of  water  is  to  transform 
an  E'  defect  and  a  nonbridging  oxygen  into  two  neutral  silanol  groups.  Simi¬ 
lar  studies  indicated  that  fibers  with  high  hydroxyl  content  were  less  prone 
to  damage  at  short  wavelengths,  although  excessive  hydroxyl  absorption  was 
observed  at  long  wavelengths  [2] .  Recently  there  have  been  many  studies  on 
increases  of  optical  attenuation  in  optical  fibers  due  to  hydrogen  permeation 
[3] .  This  loss  increase  phenomenon  is  associated  with  the  diffusion  of  hy- 
droger.  vo  the  fiber  core  and  subsequent  chemical  reactions  with  structural 
defects  in  the  glass  to  form  CH  species  and  new  defect  centers.  To  test  the 
applicability  of  hydrogen  treatment  in  defect  passivation,  we  have  compared 
the  radiation  resistance  of  a  treated  fiber  with  that  of  an  untreated  length  of 
the  same  fiber. 

In  the  first  experiment,  we  treated  a  250-m-long  section  from  a  Ge/P-co- 
doped,  graded-index  optical  fiber.  This  piece  was  placed  in  an  atmosphere 
with  a  hydrogen  partial  pressure  of  0.06  atm  at  70°C  for  seven  days.  After 
the  treatment,  dissolved  hydrogen  was  evacuated  from  the  fiber.  Figure  1 
shows  the  change  in  loss  spectrum  after  the  hydrogen  pretreatment.  The 
hydroxyl  concentration  was  found  to  have  increased  from  73  ppb  to  115  ppb. 
The  additional  loss  above  1.5  um  was  due  to  the  formation  of  POH.  Portions 
of  the  hydrogen  treated  fiber  and  an  untreated  control  fiber  were  irradiated 
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Figure  1 

Loss  spectra  of  fiber  before  and  after  hydrogen  treatment  (0.06  atm.  ?0°C. 
7  days)  showing  the  resultant  increase  in  hydroxyl  absorption. 


8 


LOSS  (dB/km) 


WAVELENGTH  (jim) 


Figure  2 

Loss  spectra  of  fibers  with  and  without  hydrogen  pretreatment  (0.06  atm, 
70®C,  7  days)  after  5  kftad  dose. 
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using  the  wCo  source  at  RADC.  The  total  dose  was  5  kRad.  Figure  2 
shows  the  loss  spectra  from  800-1600  nm  for  the  treated  and  untreated  fibers 
after  irradiation.  The  high  loss  at  shorter  wavelengths  is  expected  for  a 
phosphorus  doped  optical  fiber.  These  measurements  reveal  a  clear  difference 
between  the  treated  and  untreated  fibers  which  increases  with  increasing 
wavelength  (Figure  2).  The  attenuation  (dB/km)  is  >20%  lower  for  the  hydro* 
gen-treated  fiber  at  1600  nm.  This  modest  difference  is  not  surprising,  con¬ 
sidering  the  low  hydrogen  partial  pressure  used  in  the  pretreatment. 

In  the  second  experiment  we  treated  fibers  in  an  atmosphere  with  a  hy¬ 
drogen  partial  pressure  of  0.86  atm  at  65°C  for  seven  days  to  observe  effects 
of  higher  hydrogen  partial  pressure  and  hydrogen  posttreatment  on  radiation- 
induced  losses.  Three  samples  were  selected  for  the  hydrogen  treatment: 
both  of  the  irradiated  fibers  from  the  first  experiment  and  a  piece  of  the  as- 
drawn  fiber.  Figure  3  shows  the  loss  spectra  of  the  irradiated  fibers  after 
the  hydrogen  posttreatment.  The  data  reveal  a  dramatic  recovery  in  attenua¬ 
tion  to  10-25  dB/km  at  1100  nm  and  20-30  dB/km  at  850  nm.  They  also  re¬ 
veal  the  difference  between  the  pretreated  and  untreated  fibers  when  sub¬ 
jected  to  hydrogen  posttreatment:  the  pretreated  fiber  has  retained  its 
overall  lower  attenuation.  All  treated  fibers  and  a  control  sample  were  sub¬ 
jected  to  *5  kRad  I  irradiation.  Figure  4  compares  the  loss  spectra  after  ir¬ 
radiation  with  and  without  the  hydrogen  treatment.  We  observe  on  the  order 
of  50  dB/km  less  induced  attenuation  at  long  wavelengths  (\  >  1  um)  for  a 
hydrogen-treated  fiber.  This  is  significantly  better  than  results  from  the 
first  experiment  and  is  attributed  to  the  use  of  a  higher  hydrogen  partial 
pressure.  We  also  observed  improvement  in  radiation  hardness  as  a  result  of 
postirradiation  hydrogen  treatment.  Posttreatment  of  irradiated  fibers  by  hy¬ 
drogen  can  significantly  reduce  existing  radiation-induced  absorption  (see 
Figure  3)  and  increase  the  hardness  of  fibers  under  subsequent  irradiation. 
Figure  5  shows  the  loss  spectra  after  a  second  irradiation  for  posttreated  fi¬ 
bers.  Note  that  there  is  little  difference  between  two  posttreated  fibers  and 
the  overall  loss  is  lower  than  that  shown  in  Figure  4.  In  all  attenuation  meas¬ 
urements,  special  care  was  taken  to  avoid  the  photobleaching  effect  observed 
with  earlier  samples. 
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Figure  3 

Loss  spectra  of  fibers  retreated  (0.86  atm,  65°C,  7  days)  after  irradiation. 
Lower  curve:  with  hydrogen  pretreatment  (0.06  atm,  70°C,  7  days).  Upper 
curve:  without  hydrogen  pretreatment. 
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Figure  5 

Loss  spectra  after  a  second  4.6  kRad  dose  for  fibers  treated  with  hydrogen 
(0.86  atm,  65°C,  7  days)  following  first  irradiation.  The  upper  curve  is  the 
fiber  which  was  treated  prior  to  the  first  irradiation  (0.06  atm,  70°C, 
7  days). 
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2.2  PHOTOLUMINESCENCE  AND  RAMAN  SCATTERING 


The  radiation -induced  loss  measurements,  which  are  the  primary  charac¬ 
terization  tool  used  in  this  program,  are  generally  not  selective  enough  to 
identify  the  specific  centers  which  are  responsible  for  the  absorption.  In  or¬ 
der  to  obtain  a  more  selective  technique  for  this  purpose,  we  have  explored 
the  use  of  Raman  scattering  and  photoluminescence  in  fibers.  The  objective 
of  this  line  of  endeavor  is  to  identify  specific  color  centers  and  precursor  de¬ 
fects  and  correlate  their  concentrations  with  preform  composition  and  fiber 
processing.  This  information,  along  with  the  knowledge  of  which  centers  are 
most  deleterious  for  the  transmission  wavelengths  of  interest,  will  enable  us 
to  predict  compositions  and  processing  conditions  which  will  produce  fiber 
with  the  greatest  radiation  resistance. 

Three  phosphorus-doped  fibers  drawn  from  the  same  preform  were  inves¬ 
tigated  using  these  techniques.*  The  samples  were  short  sections  taken  from 
the  fibers  discussed  in  Section  2.1.  These  were  a  pristine  as-drawn  fiber 
(lower  curve  in  Figure  1),  a  fiber  which  was  irradiated  without  pretreatment 
(upper  curve  in  Figure  4),  and  a  fiber  which  was  hydrogen  treated  and  then 
irradiated  (lower  curve  in  Figure  4).  Raman  scattering  was  performed  on 
these  fibers  using  laser  excitation  at  wavelengths  ranging  from  482  to  676  nm. 
However,  it  was  found  to  be  relatively  insensitive  to  the  defects  of  interest, 
and  no  variations  could  be  found  among  the  fibers  or  as  a  function  of  wave¬ 
length  for  a  given  fiber.  Photoluminescence  was  found  to  be  a  very  sensitive 
probe  of  defect  centers  and  significant  differences  among  the  fibers  were  de¬ 
tected.  Figure  6  compares  the  luminescence  band  at  650  nm  for  the  three  fi¬ 
bers  before  they  were  subjected  to  short  wavelength  (X<530  nm)  excitation. 
This  emission  band  and  an  associated  absorption  band  near  630  nm  have  been 
reported  on  by  several  authors  who  attribute  it  to  drawing-  and/or  radiation  - 
induced  centers  [4].  To  make  an  absolute  comparison  of  intensities,  the 
emission  spectra  were  scaled  by  equilizing  the  intensities  of  the  Raman  scat¬ 
tering  for  the  different  fibers.  Luminescence  is  not  detectable  for  the  pris¬ 
tine  fiber  but  is  prominent  for  the  irradiated  fibers,  being  even  more  intense 
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Figure  6 

Comparison  of  luminescence  intensity  prior  to  excitation  at  482  nm  for  fibers 
subjected  to  different  treatments.  The  excitation  wavelength  was  568  nm. 
The  sharp  structure  at  left  is  Raman  scattering  which  was  used  to  scale  the 
intensities  for  the  different  fibers. 


for  the  fiber  which  had  been  hydrogen  treated  before  irradiation.  Since  the 
latter  fiber  also  exhibited  lower  losses  at  wavelengths  >1  urn  (lower  curve  in 
Figure  4),  this  inverse  relationship  to  radiation-induced  loss  implies  that  this 
defect  is  not  responsible  for  the  induced  long  wavelength  loss  in  these  fibers. 
Moreover,  it  appears  that  hydrogen  treatment  prior  to  irradiation  can  convert 
defects  which  will  result  in  long  wavelength  losses  into  the  more  benign  de¬ 
fects  which  are  transformed  into  this  luminescent  center  when  irradiated. 

The  photoluminescence  was  also  found  to  depend  upon  the  exciting  wave¬ 
length.  Figure  7  compares  the  emission  for  the  three  fibers  when  excited  at 
482  nm.  Again,  the  spectra  have  been  scaled  by  the  Raman  intensity  to  facil¬ 
itate  absolute  comparisons.  The  fibers  maintain  the  same  hierarchy  of  intensi¬ 
ties  when  excited  at  this  wavelength,  but  the  emission  is  considerably 
stronger  with  even  the  pristine  fiber  showing  intense  emission.  Additional 
features,  such  as  the  bump  at  600  nm,  have  appeared  in  the  spectra  and  may 
indicate  the  presence  of  still  other  centers.  These  have  not  been  reported  in 
the  literature  and  will  be  investigated  further. 

Prolonged  excitation  at  short  wavelengths  (X<530  nm)  produced  signifi¬ 
cant  photobleaching  in  the  irrad iatied  fibers  which  was  observed  through  vis¬ 
ual  inspection.  This  occurs  through  photoconversion  of  centers  which  is  also 
observed  through  changes  in  photoluminescence  intensity  for  all  three  fibers. 
Figure  8  illustrates  this  phenomenon  for  the  pristine  fiber.  Prior  to  short 
wavelength  exposure  there  had  been  no  significant  emission  at  650  nm  when 
excited  at  568  nm.  Afterwards  the  luminescence  band  has  become  noticeable 
although  still  much  less  intense  than  the  Raman  scattering.  Figure  9  shows 
the  same  effect  for  the  fiber  which  had  been  hydrogen  treated  prior  to  irra¬ 
diation.  The  luminescence  is  considerably  stronger,  implying  an  increase  in 
concentration  of  this  center,  although  the  fiber  loss  has  been  significantly 
reduced  by  photobleaching. 

Photoconversion  of  defects  is  a  well  known  phenomenon  in  insulators  and 
its  manifestation  as  photobleaching  has  been  extensively  studied  in  fibers. 
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Figure  7 

Comparison  of  luminescence  intensity  for  fibers  subjected  to  different 
treatments.  The  excitation  wav^'ength  was  482  nm.  The  sharp  structure  at 
left  is  Raman  scattering  which  was  used  to  scale  the  intensities  for  the 
different  fibers. 
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Figure  8 

Comparison  of  luminescence  intensity  before  and  after  482-nm  excitation  for 
pristine  fiber.  The  excitation  wavelength  was  568  nm.  The  sharp  structure 
at  left  is  Raman  scattering  which  was  used  to  scale  the  intensities. 
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Figure  9 

Comparison  of  luminescence  intensity  before  and  after  482-nm  excitation  for 
fiber  which  was  treated  with  hydrogen  prior  to  irradiation.  The  excitation 
wavelength  was  568  nm.  The  sharp  structure  at  left  is  Raman  scattering 
which  was  used  to  scale  the  intensities. 


Photoluminescence  investigations  such  as  this  one  will  be  extremely  helpful  in 
gaining  new  insight  into  types  of  defects  present  in  fibers  and  their  correla¬ 
tion  with  radiation-induced  loss  at  tong  wavelengths.  It  has  provided  addi¬ 
tional  evidence  for  our  view  that  radiation-induced  loss  can  be  reduced  by 
defect  conversion  (passivation)  before  as  well  as  after  irradiation.  The  cen¬ 
ter  responsible  for  the  650-nm  luminescence  band  is  an  example  of  a  relatively 
benign  defect  which  can  be  formed  by  hydrogen  treatment  from  more  deleteri¬ 
ous  defects.  This  investigation  will  be  extended  to  analyze  the  effects  of 
other  compositions  and  processing  parameters. 


3.  DEFECT  MINIMIZATION 


3.1  MOLECULAR  DYNAMICS  SIMULATION 

The  molecular  dynamics  simulation  uses  computer-generated  glasses  to 
investigate  the  types  and  abundances  of  structural  defects  present  in  glass 
compositions  used  for  fibers.  The  objective  is  to  determine  the  effect  on  the 
defect  structure  of  silica  glass  of  conventional  dopants  as  well  as  dopants 
which  serve  as  recombination  centers  for  the  electrons  and  holes  produced  by 
radiation.  By  gaining  a  fundamental  understanding  of  what  types  and  con¬ 
centrations  of  defects  are  formed  by  dopants,  we  hope  to  predict  which 
dopants  or  combination  of  dopants  will  result  in  the  greatest  radiation  hard¬ 
ness  while  providing  the  required  refractive  index  modification  or  recombina¬ 
tion  centers. 

Molecular  dynamics  generates  the  glass  structure  by  calculating  the 
pairwise  interionic  forces  of  a  large  array  of  ions  and  uses  these  to  determine 
the  velocity  and  position  of  the  ions  at  a  later  time.  The  new  positions  are 
in  turn  used  to  recalculate  the  interionic  forces-  This  cyclical  process  is  re¬ 
peated  at  short  time  intervals  at  a  high  temperature  until  the  structure  stabi¬ 
lizes.  The  glass  is  then  cooled  to  room  temperature  by  removing  kinetic  en¬ 
ergy  and  suitable  analysis  programs  are  applied  to  examine  the  local  structure 
of  this  computer-generated  glass.  Analysis  programs  currently  available  cal¬ 
culate  the  radial  distribution  function  of  each  type  of  ion,  the  average  diffu¬ 
sion  rate,  the  velocity  autocorrelation  function,  the  distribution  of  neighbors, 
the  fraction  of  bonds  rearranged  or  unchanged  as  a  function  of  time,  and  the 
time-dependent  radial  distribution  fraction.  The  program  used  to  generate 
the  glass  requires  input  parameters  which  describe  the  interionic  forces  of  all 


the  species.  In  particular,  repulsion  parameters  are  required  for  all 
components;  if  they  are  not  available  in  the  literature  these  parameters  must 
be  chosen  empirically  such  that  interatomic  distances  agree  with  average  val¬ 
ues  determined  experimentally. 

The  computer  program  that  generates  the  glass  structures  required  par¬ 
tial  rewriting  because  (1)  the  direct  access  files  used  by  the  program  as  it 
had  been  written  and  initially  used  were  too  large  to  retrieve  from  the  ar¬ 
chives  with  the  current  disk  usage,  and  (2)  the  FORTRAN  compiler  that  had 
been  used  previously  had  some  serious  problems  and  was  replaced  by  a  more 
universally  compatible  one.  With  the  assistance  from  a  systems  programmer 
from  the  GTE  Technical  Computation  Center,  the  program  was  modified  to  be 
compatible  with  the  new  compiler  and  use  smaller  files  as  well  as  run  more  ef¬ 
ficiently  . 

The  first  application  of  the  simulation  program  was  to  pure  Si02  which 

was  modeled  as  a  standard  to  which  more  complex  systems  can  be  compared. 

The  glass  was  initially  equilibrated  at  6000  K  for  12-16  picoseconds  (ps)  and 

then  cooled  gradually  to  room  temperature  in  another  16  ps.  The  time  scale 

-3 

for  each  step  in  the  simulation  is  10  ps  and  each  step  requires  =17  sec  of 
CPU  time  on  the  IBM  3081.  Another  computer  program  was  written  to  analyze 
the  structure  of  the  simulated  glass  by  calculating  the  radial  distribution 
function  and  the  average  coordination  number  as  a  function  of  distance.  This 
program  was  applied  to  SiOj  and  the  results  are  illustrated  in  Figure  10. 

Studies  were  made  to  ensure  that  the  glass  was  being  equilibrated  at 
6000  K  for  a  time  sufficiently  long  for  the  structure  to  converge.  This  was 
done  by  holding  the  glass  for  longer  periods  of  time  at  the  high  temperature 
before  cooling  to  room  temperature.  The  radial  distribution  function  was  then 
calculated  and  examined  for  change  as  a  function  of  additional  time  at  high 
temperature.  This  investigation  revealed  that  about  16  ps  of  simulation  at 
6000  K  followed  by  16  ps  of  cooling  to  room  temperature  are  needed  to  reach 
an  equilibrated  structure,  i.e.  one  that  does  not  change  significantly  with 
further  time  at  6000  K. 


. — r  |V 

l.ZSS&.J.p jw 
!  V  // 


OJ  0.4 
DISTANCE 


RAM  At.  OtSTWaunOM  RWCTIOfl  ■■  SI-0 

Mft IRAQ*  COORDINATION  NUMMR - SI-0 


— - St-Si 

041  — - $141  ■ 


Figure  10 

Radial  distribution  function  and  average  coordination  number  as  a  function  of 
distance  for  fused  silica. 
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Based  upon  these  results  for  SiO.,,  pure  Ge00  glasses  and  GeO-, -c  ped 
silica  glasses  were  generated  by  starting  with  an  equilibrated  SiC^  structure, 
adjusting  the  density  and  therefore  the  box  size,  and  allowing  the  structure 
to  equilibrate  for  an  additional  20  ps  before  cooling.  For  the  Ge-containing 
glasses,  the  same  parameters  for  oxygen  were  used  as  in  the  pure  SiOj 
glass.  The  repulsive  parameters  for  Ge  were  found  by  adjusting  their  values 
until  the  experimental  interionic  distances  were  obtained.  The  pure  Ge07 
glass  was  simulated  first  in  order  to  achieve  this.  The  Ge-0  distance  in  the 
simulated  pure  GeO?  glass  is  0.174  nm.  the  same  as  in  x-ray  diffraction  stud¬ 
ies.  The  Ge-Ge  and  0-0  distances  deviated  by  0.01  to  0.02  nm  from  the  ex¬ 
pected  values.  These  could  not  be  adjusted  without  drastically  changing  the 
values  for  oxygen.  This  problem  in  pure  Ge02  may  reflect  the  covalent  na¬ 
ture  of  the  Ge-0  bond  and  illustrates  the  shortcomings  of  using  a  purely 
ionic  potential  to  generate  the  glass  structure.  However,  in  glasses  contain¬ 
ing  primarily  SiOj,  the  Ge-Ge  interactions  should  be  minimal  and  this  discre¬ 
pancy  is  not  expected  to  have  a  significant  effect  upon  the  structure. 

As  the  model  system  for  Ge-doped  silica,  a  glass  with  the  composition 
10  mole*  GeOj  -  90  mole*  SiC^  was  generated.  The  simulation  was  done  by 
starting  with  the  equilibrated  SiC^  structure  and  replacing  6  of  the  60  Si  ions 
by  Ge  ions.  The  Ge  parameters  used  were  those  obtained  in  the  above  simu¬ 
lation  of  pure  GeO  . .  The  structure  of  the  binary  glass  was  then  equilibrated 
at  6000  K  for  20  ps  before  cooling  to  room  temperature.  To  obtain  a  statisti¬ 
cally  valid  sample,  10  cooling  runs  were  done.  Each  run  started  with  the  last 
configuration  of  the  previous  6000  K  run  and  was  allowed  to  continue  for  an 
additional  2000  steps  (2  ps)  at  6000  K  before  cooling.  After  cooling  to  300  K 
and  equilibrating  for  1000  steps,  every  tenth  configuration  was  saved  so  that 
from  the  ten  runs,  1000  configurations  could  be  examined.  The  radial  distri¬ 
bution  function,  the  average  coordination  number,  and  the  probability  that  a 
Si  ion  has  3,4,  or  5  oxygen  ions  within  a  certain  distance  were  determined. 
Furthermore,  differences  between  'average'  Si  ions  and  Si  ions  near  Ge 
ions”  could  be  distinguished. 
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Figure  11 

Radial  distribution  function  and  average  coordination  number  as  a  function  of 
distance  for  the  10  mole°o  GeO^ -doped  silica  glass.  The  vertical  scale  has 
been  expanded  and  the  peak  of  the  rdf  curve  lies  off  scale 
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Figure  12 

The  average  coordination  number  as  a  function  of  distance  for  the  10  mole'o 
GeOj-doped  silica  glass.  The  ordinate  has  been  expanded  over  that  of 
Figure  It  to  illustrate  the  differences  at  distances  beyond  the  median  Si-0 
separation . 


Examination  of  the  coordination  number  at  distances  greater  than  the 
median  Si-0  distance  but  within  the  first  coordination  sphere  shows  that  a 
number  of  Si  ions  have  an  oxygen  coordination  number  greater  than  four.  As 
seen  in  Figure  11,  at  distances  less  than  0.215  nm  the  oxygen  coordination 
around  Si  ions  near  Ge  ions  is  greater  than  that  of  "average"  Si  ions.  At 
distances  greater  than  0.215  nm,  the  "average"  Si  ions  have  a  larger  coordi¬ 
nation  number.  This  can  be  seen  clearly  in  Figure  12  which  shows  the  Si-0 
coordination  number  as  a  function  of  distance  on  an  expanded  scale.  The 
same  effect  is.  reflected  in  the  radial  distribution  function  (RDF).  In  Fig¬ 
ure  11  only  the  edges  of  the  main  RDF  peaks  are  shown,  the  peaks  lie  ■  off- 
scale  vertically.  The  nearest  neighbor  Si-0  RDF  for  Si  ions  near  Ge  ions 
shows  a  "tighter"  structure  than  for  "average"  Si  ions.  At  distances  greater 
than  0.20  nm  from  Si  ions,  the  probability  of  finding  an  oxygen  ion  is  greater 
for  'average"  Si  ions  than  for  Si  ions  near  Ge  ions.  These  would  be  the 
fifth  oxygen  surrounding  a  central  Si  ion.  Thus  the  Si04  tetrahedra  occa¬ 
sionally  have  an  additional  oxygen  ion  crowding  them.  Another  way  of  look¬ 
ing  at  these  results  is  that  the  coordination  number  of  Si  ions  near  Ge  ions  is 

closer  to  four  than  that  of  "average"  Si  ions  because  of  the  narrower  distri¬ 

bution  in  distances  to  oxygens.  This  is  also  represented  in  the  probability 
that  Si  ions  will  have  an  oxygen  coordination  number  of  five  as  shown  in  Ta¬ 
ble  1.  The  effect  is  particularly  noticeable  at  distances  of  0.215-0.24  nm 
from  Si  ions. 

Thus  far  this  analysis  indicates  that  defects,  in  terms  of  stretched  Si-0 
bonds  near  Ge,  are  not  observed  in  this  glass.  Defects  produced  by  Ge  dop¬ 
ing  which  are  susceptible  to  radiation  must  be  of  another  type.  They  could 
be  found  in  Ge-0  bonds  themselves  or  in  "average"  Si-0  bonds  which  differ 
from  those  in  pure  S^.  From  x-ray  studies  of  pure  GeOj,  it  has  been 
found  that  there  is  a  very  narrow  distribution  of  Ge-O-Ge  intertetrahedral 
angles  and  the  randomness  in  the  glass  structure  is  attributed  to  tetrahedra 

twisting.  Whether  this  can  be  extended  to  Si02*Ge02  glass  has  yet  to  be  de¬ 

termined.  If  this  is  the  case,  then  constraints  may  be  placed  on  the  network 
which  are  not  seen  in  pure  Si02-  The  defect  structure  would  then  be  more 


TABLE  1 


PROBABILITY  OF  Si  HAVING  N  OXYGEN 
NEIGHBORS  AS  A  FUNCTION  OF 
NEIGHBOR  DISTANCE  D 


A.  All  Si  Ions 


PROBABILITY 

0.146  <  D  <  0.186 
(nm) 

0.146  <  D  <  0.215 
(nm) 

0.146  <  D  <  0.240 
(nm) 

0.024 

0.973 

0.003 


0 

0.952 

0.048 


B.  Si  Ions  With  a  Ge  Neighbor  at  Distance  Within  Range  0.25>0.4  nm 


_  PROBABILITY  _ 

N  0.146  <  D  <  0.186  0.146  <  D  <  0.215  0.146  <  D  <  0.240 
(nm)  (nm)  (nm) 


.939 

.061 


complicated  than  just  abnormal  ionic  distances  and  might  include  ring 
structures. 

Oetailed  structural  analysis  of  this  type  will  also  be  performed  on  pure 
SiOj  and  Ce02  glasses.  Comparison  of  those  results  to  the  ones  presented 
here  for  the  doped  glass  will  assist  us  in  determining  the  types  of  defects 
generated  when  silica  and  germama  are  mixed.  This  analysis  is  in  progress. 
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To  continue  the  molecular  dynamics  study  of  optical  fiber  glasses,  a 
SiOj-GeOj'F  glass  is  being  simulated.  The  repulsive  parameter  being  used 
for  fluorine  is  that  of  Tosi  and  Fumi  for  fluoride  salts  [5] .  The  initial  struc¬ 
ture  is  that  of  an  equilibrated  SiOj-GeOj  glass  with  the  appropriate  substitu¬ 
tions,  and  the  structure  is  being  allowed  to  equilibrate  for  20  ps  at  6000  K 
before  cooling. 


3.2  DRAWING  INDUCED  DEFECTS 

In  this  aspect  of  the  program  the  influence  of  drawing  conditions  on  the 
radiation  sensitivity  of  fiber  is  investigated.  Four  fibers  were  drawn  from 
the  same  preform  at  extremes  of  the  interdependent  parameters  of  tempera¬ 
ture,  speed,  and  tension.  Examination  of  the  loss  spectra  before  and  after 
irradiation  will  indicate  the  role  of  the  fiber  drawing  process  in  generating 
defects  which  become  color  centers. 

A  graded  index  preform  was  fabricated  for  this  investigation  of  drawing- 
induced  defects.  The  preform  had  a  constant  fluorine  doping  level  across 
both  clad  and  core  and  a  Ge  doped  core.  Four  800-m  sections  of  optical  fiber 
were  drawn  under  the  four  different  sets  of  drawing  conditions  listed  in  Ta¬ 
ble  2.  The  furnace  temperature  and  draw  speed  were  selected  to  intentionally 
vary  the  induced  defect  concentration  among  the  fibers.  Figures  13-16  illus¬ 
trate  the  loss  spectra  of  these  fibers  measured  from  600-1600  nm  before  and 
after  a  T  radiation  dose  of  »20  kRad.  The  oscillatory  structures  in  Figures 
14  and  16  are  probably  due  to  weak  guiding  resulting  from  the  depressed  in¬ 
ner  clad.  We  intend  to  study  in  greater  detail  the  differences  in  fiber  atten¬ 
uation  prior  to  irradiation  and  their  relationship  to  drawing-induced  defects 
using  additional  techniques  such  as  Raman  scattering  and  photoluminescence 

Table  2  also  summarizes  the  loss  results  at  850  nm  before  and  after  irra¬ 
diation.  The  fibers  drawn  at  the  lower  furnace  temperature  showed  an  over- 
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FIBER  B 


figure  14 

Los*  spectrum  of  fiber  B  before  and  after  22  kRad  irradiation 
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Figure  15 

Loss  spectrum  of  fiber  C  before  and  after  22  kRad  irradiation. 
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Figure  16 

Loss  spectrum  of  fiber  D  before  and  after  24  kRad  irradiation. 
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TABLE  2 


SUMMARY  OF  RESULTS  FOR  THE 
INVESTIGATION  OF  THE  EFFECT  OF  DRAWING 
PARAMETERS  ON  RADIATION-INDUCED  LOSS 


all  lower  radiation-induced  loss,  especially  at  longer  wavelengths  (si. 3  um) . 
The  loss  values  at  850  nm  agreed  with  in  sJtu  results  obtained  at  RADC 
shown  in  Figure  17.  It  is  generally  believed  that  fiber  drawn  at  a  lower 
temperature  has  a  lower  fictive  temperature  and,  consequently,  a  lower  defect 
concentration.  This  is  in  agreement  with  our  experimental  results.  The  ef¬ 
fect  of  draw  speed  is  less  clear  and  further  study  will  be  necessary  in  order 
to  understand  it.  The  results  of  this  investigation,  specifically  the  advan¬ 
tages  of  fluorine  doping  in  the  core  and  the  use  of  lower  draw  temperatures, 
are  also  being  applied  to  subsequent  work  in  the  other  approaches. 


3.3  STRESS- 1 NOUCED  DEFECTS 


Thermomechanical  stresses  arise  in  optical  fibers  because  of  the  different 
thermal  characteristics  of  the  glasses  used  for  the  core  and  the  cladding.  In 
general,  this  produces  a  tensile  stress  in  the  more  heavily  doped  core  causing 
atomic  bonds  to  lengthen  and  possibly  rupture,  thereby  creating  defects. 
This  study  is  intended  to  assess  the  impact  of  thermomechanical  stress  on  ra¬ 
diation  hardness  and  identify  ways  to  reduce  the  number  of  stress-induced 
defects . 

To  examine  the  effect  of  core/clad  interfacial  stress  on  defect  concentra¬ 
tions,  we  prepared  two  preforms  with  different  clad  compositions  which  also 
kept  the  index  profile  the  same.  The  first  one,  having  a  pure  silica  clad  and 
Ge02/F  codoped  core,  provides  a  high  core/clad  interfacial  stress  due  to  dif¬ 
ferences  in  thermal  expansion  and  glass  transition  temperatures.  The  second 
one,  having  a  GeO^/F  codoped  clad  and  GeOj/F  codoped  graded-index  core, 
minimizes  the  core/clad  stress.  In  both  preforms  the  fluorine  doping  was 
maintained  at  a  constant  level  and  the  index  of  the  clad  was  matched  to  that 
of  the  substrate  tube.  Electron  microprobe  results  verify  that  the  desired 
structure  was  achieved  and  reveal  high  concentration  gradients  for  both  Ge 
and  F  across  the  core/clad  interface  in  the  first  preform,  in  contrast  to 
slowly  varying  profiles  for  the  second  preform.  We  have  drawn  both  pre¬ 
forms  into  fibers  using  optimized  draw  parameters  as  determined  by  the  re¬ 
sults  of  the  above  investigation  of  drawing-induced  defects  (Section  3.2). 
Two  draw  speeds  were  used  for  each  preform. 

Table  3  summarizes  results  for  the  investigation  of  the  effect  of  stress 
reduction  on  radiation-induced  loss.  Fibers  designated  HA  and  HB  were 
drawn  from  the  preform  with  a  high  core/clad  interfacial  stress,  while  fibers 
LA  and  LB  were  from  a  low  stress  preform.  All  these  fibers  show  excellent 
radiation  resistance.  The  induced  losses  at  both  850  nm  and  1300  nm  are 
quite  low  compared  to  fibers  with  only  GeC^  doped  cores.  Figures  18  to  21 
compare  loss  spectra  of  these  fibers  before  and  after  a  dose  of  =20  kRad.  To 


TABLE  3 


SUMMARY  OF  RESULTS  FOR  THE 
INVESTIGATION  OF  THE  EFFECT  OF  STRESS 
REDUCTION  ON  RADIATION-INDUCED  LOSS 


Fiber 

Designation 

Draw 

Speed 

(m/s) 

Loss  a 

850  nm  (dfl 

/km) 

Loss  at  1300  nm  (d 

B/km) 

Before 

Irradiation 

After 

irradiation 

Induced 

Loss 

Before 

Irradiation 

After 

Irradiation 

Induced 

Loss 

HA 

wm 

2.43 

34.2 

31J 

0.80 

1.46 

0.66 

HB 

Q 

2.63 

34.0 

31.4 

0.92 

3.84 

2.92 

LA 

IjjgR 

2.42 

33.2 

30.8 

0.84 

3.14 

2.30 

LB 

■9 

2.58 

29.9 

m 

0.92 

2.38 

1.46 

further  clarify  the  stress  effect,  we  used  the  differential  modal  attenuation 
(DMA)  technique  to  monitor  the  attenuation  of  a  specific  group  of  modes. 
Figures  22  and  23  show  DMA  results  at  1300  nm  for  fibers  drawn  from  pre¬ 
forms  having  high  and  low  stress,  respectively.  In  general,  DMA  data  dem¬ 
onstrate  higher  losses  for  higher  modes.  The  DMA  measurements  do  not  cor¬ 
relate  well  with  the  attenuation  measurements,  however,  since  mode  filters 
were  used  for  the  latter.  Mode  filters  remove  the  higher  modes,  and  thus 
the  loss  results  only  reflect  attenuation  of  lower  modes.  Before  irradiation 
the  stress  causes  a  steady  increase  in  attenuation  with  increasing  mode  num¬ 
ber  (Figure  22).  For  a  low  stress  fiber  (Figure  23)  the  increase  is  slight 
and  gradual.  After  irradiation,  we  do  not  observe  any  significant  difference 
between  fibers  with  different  stresses  in  attenuation  and  DMA  data. 
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Figure  19 

Loss  spectrum  of  high  stress  fiber  HB  before  and  after  21.1  kRad  irradiation. 


Figure  22 

DMA  results  before  and  after  irradiation  of  fibers  with  high  stress. 


Figure  23 

OMA  results  before  and  after  irradiation  of  fibers  with  low  stress. 
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Prior  to  this  investigation,  there  was  no  known  information  on  the  spatial 
distribution  of  radiation-sensitive  defects  produced  by  core/clad  interfacial 
stress.  These  results  suggest  that  these  stresses  introduce  defects  only  near 
the  interface,  resulting  in  higher  attenuation  for  higher  modes.  Since  the 
higher  modes  do  not  propagate  a  great  distance  in  a  fiber,  however,  stress 
reduction  does  not  contribute  materially  to  radiation  hardness.  Moreover,  the 
absence  of  significant  differences  in  the  attenuation  of  the  higher  modes  fol¬ 
lowing  irradiation  implies  that  the  defects  produced  by  stress  either  are  not 
great  in  number  or  are  not  very  sensitive  to  I  radiation. 


3.4  DOPANT  REDUCTION 

Dopants  influence  the  concentration  of  defects  in  as-drawn  fibers,  as 
they  do  in  bulk  glass,  by  altering  the  glass  structure.  To  evaluate  the  ef¬ 
fect  of  dopant  reduction  on  radiation  resistance,  we  have  fabricated  three 
preforms  with  different  Ge  and  F  dopant  concentrations  in  the  core.  These 
preforms  have  concentration  ratios  of  1/0.86/0.73  for  germanium  and  1/0. 5/0 
for  fluorine.  The  reduction  was  introduced  in  such  a  way  as  to  preserve  the 
index  profile  and  minimize  the  core/clad  stress.  Fibers  drawn  from  this  se¬ 
ries  of  preforms  will  be  irradiated  at  the  RADC  I  source  and  fully  character- 
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